Localized Rayleigh Instability in Evaporation Fronts 
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A qualitatively different manifestation of the Rayleigh instabihty is demonstrated, where, instead 
of the usual extended undulations and breakup of the hquid into many droplets, the instability is 
localized, leading to an isolated narrowing of the liquid filament. The localized instability, caused 
by a nonuniform curvature of the liquid domain, plays a key role in the evaporation of thin liquid 
films off solid surfaces. 
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The Rayleigh instability of slender liquid bodies, 
driven by surface tension, is part of everyday experi- 
ence and has been systematically studied for well over 
a century [H, @|. Its manifestations are diverse, ranging 
from the breakup into droplets of inviscid, viscous, and 
viscoelastic liquidjets and bridges [H to the pearling of 
fluid membranes [3| . The liquid bodies considered in ear- 
lier studies were usually translation-invariant along their 
long dimension, resulting in unstable modes which were 
extended [l], Q or had a steadily propagating front [1] . 
The propagation of forced perturbations from a fixed noz- 
zle were studied as well jlj]. In this Letter we investigate 
a qualitatively different scenario of the Rayleigh insta- 
bility, in which the translation invariance is broken by a 
nonuniform curvature, and show that the fastest-growing 
mode of this instability is localized. The localization is 
reminiscent of the problem of a quantum particle moving 
inside a curved stripe '^% ifj] . When the latter system is 
transformed into an effective one-dimensional problem, 
an attractive potential emerges, whose minimum is lo- 
cated at the point of maximum curvature, giving rise to 
bound (localized) states. 

Thin liquid bodies are abundant in phenomena related 
to wetting of solid substrates [1, • Two processes are 
particularly relevant to the current work: (i) the dewet- 
ting of a nonvolatile, nonwetting film 0,111; and (ii) 
the evaporation of a volatile, totally wetting film fiol- 13|. 
Both processes exhibit the kinetics of a first-order tran- 
sition, where dry domains [in (i)], or domains covered by 
a molecularly thin liquid [in (ii)], nucleate and grow into 
a much thicker film. Importantly, in both processes the 
dewetting front has a long, slender rim of excess fluid Q . 
The growing domains may have a stable circular bound- 
ary [8|, llC| or evolve through elaborate instabilities and 
patterns [10|-tlJ|. An analogy has recently been drawn 
between the pattern formation in the volatile case and 
Saffman- Taylor viscous fingering in a Hele-Shaw cell hav- 
ing a time- varying thickness |l5l | . The additional dynam- 
ics of the liquid rim at the domain boundary, however, 
crucially affects the selection of patterns in experiments, 
e.g., the doublon pattern shown in Fig. [Tl which is un- 
common in Saffman- Taylor fingering ll|, Il5| . 



We model the droplet 's rim as a long, curved liquid 
strip of uniform width w. It is parametrized using the 
locally orthogonal triad (s, u, z) as depicted in Fig. [21 
where u e [~w/2,w/2] and z e [0,h], h{s,u) being the 
local liquid height. The rim is in contact with the wet 
and dry domains at m = —w/2 and 'w/2, respectively. 
The local curvature of the centerline (s, 0, 0) is denoted by 
k{s) and taken as positive when the rim curves away from 
the dry domain. (In Fig.[5]K < 0.) We assume separation 
of time scales between the growth of the finger and the 
faster development of the rim instability, as observed in 
the evaporation experiments. Thus, although the rim 
itself is a dynamic effect, w and k(s) can be taken as 
time- independent . 

Within the lubrication approximation (|V/i| ^ 1) flG} . 
the equation of motion for the liquid is 



3r] 



(1) 



where p is the liquid pressure, ij its viscosity, and spa- 




FIG. 1: Doublon patterns observed in the evaporation of wa- 
ter off a clean mica surface • Fingers of a molecularly thin 
water film grow into a thicker film. A local indentation at 
the finger tip leads to splitting and the formation of a liquid 
spine that subsequently breaks into droplets. The fingers are 
hundreds yum wide. 
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FIG. 2: Schematic view of the system and its parametriza- 
tion. A curved liquid rim of uniform width w, nonuniform 
curvature k{s), and height profile z = h{s,u), lies at the in- 
terface between a thick liquid film (wet domain) and a much 
thinner one (dry domain). 

tial derivatives are in the (s, u) plane. In p we include 
contributions from the Laplace pressure and disjoining 
pressure n(/i) (i.e., surface interactions) ^3], 

p^-jV^h~U{h), (2) 

where 7 is the surface tension of the liquid. Substituting 
in Eqs. ^ and ([2]) a small perturbation of the steady 
profile, h = ho+ijj, linearizing in -0, and neglecting spatial 
derivatives of /iq, we obtain 

5t0 = OV, O - -ro^o"^ (V^ + 2g2v2) . (3) 

In Eq. ([3]) the Laplacian is given by ^ 

<7(s,u) = [l + uK(s)]^ (4) 

and the following length and time scales appear: 

= [2j/U'{ho)]'/'. r„-i = 12r^^/ihl[U'iho)r). (5) 

In the systems under consideration the disjoining pres- 
sure is usually governed by van der Waals interactions 
0, whereby W = H/{2T:h^), H - lO-^^ erg being the 
Hamaker constant. The length b = {H/^y^'^ is invari- 
ably of molecular scale, ~ 1 nm. Hence, the value of 
^ h^/b is primarily determined by the rim thickness, 



ho, which has very different values in the dewetting and 
evaporation processes. In the dewetting case of Ref. ^ 
ho is about 50 iim, leading to q^^ of order meters and 
an unphysically long T'^^ . In the evaporation process of 
Ref. J^], by contrast, ho is of order 10 nm, yielding q^^ 
of micron scale and Fq ^ of order 10^^ s (though the lat- 
ter is highly sensitive to the thickness). The rim width, 
w, is of mm scale in both cases; thus, wqo <C 1 for the 
dewetting case, and wqo ^ 1 for the evaporation one. 
As we shall presently see, instability of the rim requires 
wqo ^ 1, which clarifies the strikingly different dynamics 
observed in the two processes. 

A key feature of the system is that the rim separates 
domains of differing properties (Fig. This implies 
asymmetric boundary conditions at u = ±w/2 and, con- 
sequently, sensitivity of the results to the curvature direc- 
tion (sign of k). At the boundary with the thick liquid 
film we impose a fixed height and a vanishing surface 
curvature, 

u = ~w/2 : ip = 0, VV = 0, (6) 

thus ensuring that the pressure [Eq. ([^] changes continu- 
ously between the rim and the film [18| . At the boundary 
with the dry domain we assume a fixed contact angle and 
a vanishing outward current, 

u^w/2: duip = 0, 9„VV = 0. (7) 

The operator O, as defined by Eq. ([3]) and the bound- 
ary conditions ([6]) and (O , is hermitian. In general it does 
not commute with the Laplacian because of the spatially 
varying g{s,u). In cases where O and do commute, 
the spectrum of O (denoted by F) can be written in terms 
of that of (A) as F = -Fog(7^A(A + 2gg), and is thus 
bounded from above by Fq. Hence, in such cases the 
fastest growing mode has (at the most) a rate Fq as given 
by Eq. (P. 

In the simple case of a straight rim [l3j, k = 0, we 
have [O, V^] = 0. The eigenmodes of O in this case are 
extended, 

= Ae'«^+^'-*sin[fc„(m-u;/2)], (8) 

where A is an arbitrary amplitude, q the wavenumber 
along the rim, fc„ — 7r(n — l/2)/w {n = 1,2,...) the 
wavenumber in the transverse direction, and 

r„, - Fo(g2 + kl)[2ql - {q^ + kl)]/qt (9) 

is the growth rate. Unstable modes, having Tnq > 0, are 
obtained for a sufficiently wide rim, wqo > 7r/(2-\/2); get- 
ting fastest growing modes of finite wavelength requires 
the slightly stricter condition wqo > 7r/2. These modes 
have wavenumbers qn = {qo ~ ^n)^^^ ^ind the maximum 
growth rate F = Fq. In the case of a circular rim of 
fixed curvature, k = kq ^ 0, despite the nonuniform 
metric g = {\ + kqu)^, O and still commute due to 
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rotational symmetry. Hence, a finite uniform curvature 
cannot accelerate the instability beyond the rate Tq. 

We now address the interesting and practically rele- 
vant case of a nonuniform curvature k{s). Motivated by 
the evaporation experiments (Fig. [T] and [lol - [l3| ) and in- 
spired by the quantum bound states (i.e., localized eigen- 
modes of V^) found in a similar curved geometry Q, we 
look for a localized unstable mode of O, whose growth 
rate exceeds Fq. To this end we employ a variational ap- 
proach, which sets a lower bound for the maximum rate 
according to 



Maximizing F with respect to a yields 



> f [v^] = (V'|0|V)/(V'IV'), 



(10) 



■0(5, u) being any trial function that satisfies the bound- 
ary conditions ([6|) and d?]). 

For a general form of k{s) it is difhcult to construct 
a good variational wavefunction which will satisfy the 
boundary conditions. To simplify the analysis we assume 
that the curvature is both small and slowly varying in 
space. In this regime the wavefunction may be locally 
approximated by that of a circular rim having curvature 
k{s). This observation leads to the following choice of 
variational wavefunction: 



-ipisju) — A(/)(s, u)e* 



2 / 2 

/a 



where 



w w 2 
4>{s,u) = {1 - k(s) — [fci(M + — ) ]}sin[A;i(u- 
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(11) 
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w , 



w ^ ^ 

-[2K{s) — ki{u + -)] cos[k^{u + -)]. (12) 

TT Z Z 

Here ki = tt/{2w), Qi = {qq — kiy/^, and the localiza- 
tion length a serves as a variational parameter 1^ 20| . 
Equation has been obtained from the asymptotic 
form of the eigenmodes for a circular rim in the limit of 
small curvature, while replacing kq with the slowly vary- 
ing curvature, k{s). For the sake of concreteness let us 
take 



(13) 



Our calculation is performed to leading order in two small 
parameters: 6 = KqW and e — {qia)~^. (Note that the 
overall turn of the rim, KqCt, may still be appreciable.) 
The smallness of e, implying that the localized mode ex- 
tends many wavelengths away from the tip, is an ansatz 
to be confirmed below. In addition, we assume qow 3> 1 
(and thus qiw ^ 1) to be safely inside the unstable 
regime. The trial function of Eqs. (fTT|) and (|12p satis- 
fies the boundary conditions for ip{—w/2) and duip{w/2) 
exactly, while those for W'^ip{—w/2) and 9„V^'0(^/2) are 
violated only at the orders and Se^, respectively. 
Within this approximation we obtain from Eq. (|10p 



f (a) = Fo 
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(14) 



= ~\/2w I {kqg 
> f (a*) = Fo ( 1 
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(15) 



Thus, provided that the rim is negatively curved {kq < 
to get a* > 0), the fastest-growing mode is localized and 
achieves a growth rate larger than Fq. The required cur- 
vature direction is in accord with the experiment, where 
the instability occurs at the dry fingertips ([nj and Fig. 
[1]). Note that e = {qia)^^ ^ \Ko\a/{qiw) <C 1, and 
Fmax/Fo — 1 ~ ^ 1, which is consistent with the 
aforementioned approximations. The main result, Eq. 
(Uni), seems robust to the choice of trial function [20|. 

Two issues have remained unspecified in the discussion 
above. First, for our choice of k{s) O is symmetric in ±s 
and, therefore, the actual unstable mode must have a 
definite parity [i.e., the e*'^^'' factor in Eq. ((TT|) should be 
replaced by either cos{qis) or sin(gis)]. In the even case 
the perturbation is maximum at the tip, whereas in the 
odd case it has a node there. Within our approximation 
the difference in F between the two functions is minute. If 
the overall turn is relatively small, \Koa\ < S^^^, we have 
a > (7, and the even perturbation is found to be slightly 
faster, with (f even -f odd )/Fo = [8ql4ay{qlw'')]e-''"''\ 
An example of an even unstable mode is shown in Fig. 
[31 In the opposite case of \Koa\ > S^^"^ we have a < a, 
and the growth rate of the odd perturbation is slightly 
higher, (fodd - fevcn)/Fo - 2iqi/qo)^e-i"^"/\ 

The second question concerns the sign of the localized 
instability — whether it increases the liquid height at 
s = toward the formation of a droplet, or decreases it 
toward pinch-off. A rigorous answer requires nonlinear 
analysis that lies beyond the scope of the current work. 
Yet, since the entire liquid film is unstable against evapo- 
ration, it is plausible to expect that the rim should shrink 
at its tip. This is also the direction required to account 
for the experimentally observed doublon patterns (Fig. 
[1]). By Darcy's law the velocity of the interface increases 
with the local film height. Hence, a local indentation 
at the tip of a dry finger advances more slowly than its 
shoulders, leading to the shape of a finger split in two by 
a narrow liquid spine. 

This work has two direct experimental implications. 
The first relates to the stable radial growth of domains, 
found in the dewetting of nonvolatile liquid films 
vs. the unstable pattern formation observed in the evap- 
oration of volatile ones [lol - [l3j | . We have shown that the 
qualitatively different dynamics in these two types of ex- 
periment can be related to the stability vs. instability of 
the accumulated liquid rim at the domain boundary. The 
second implication concerns the mechanism behind the 
patterns selected in the evaporation process. The newly 
demonstrated effect — a localized Rayleigh instability 
driven by the surface tension of a nonuniformly curved 



4 



o 





FIG. 3: Localized instability of a curved liquid rim. (a) Longi- 
tudinal profile of the height perturbation along the centerline, 
■i/)(s,0). (b) Transverse profile at the tip, ipiQyu). (c) (color 
online). Two-dimensional topography, il){s,u). The parame- 
ters used are hq/w = —0.02, a/w = 7, and qow = 2. 



liquid domain — is essential for the pattern formation as 
it suppresses the SafFman- Taylor instability at the finger 
tips. More broadly, since the relation between inhomo- 
geneity and localization is far more general (encountered, 
e.g., in the effect of defects on electron states in a solid), 
we expect related localization effects to emerge in other 
scenarios of the Rayleigh instability. Indicating these 
scenarios calls for further investigation. Such localized 



instabilities, for example, may offer new possibilities to 
control via curvature the precise location where a liquid 
filament is to agglomerate into a drop or pinch off. 
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The picture used in Fig. [1] is courtesy of S. G. Lipson. 
The boundary condition ([6]) implies that, unlike other 
Rayleigh instabilities, the rim can exchange liquid with 
the film and its mass is not conserved. 
One can treat qi as another variational parameter, yet 
this does not significantly improve the results. 
Identical results are obtained using a different trial func- 
tion, (j> — sin[A;i(u + w/2)] times a 3rd-order polynomial 
in u, whose coefficients are so chosen as to satisfy the 
boundary conditions Q and ((TJ for fixed k. 



